Anomalous scaling and super-roughness in the growth of CdTe polycrystalline films 
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CdTe films grown on glass substrates covered by fluorine doped tin oxide by Hot Wall Epitaxy 
(HWE) were studied through the interface dynamical scaling theory. Direct measures of the dy- 
namical exponent revealed an intrinsically anomalous scaling characterized by a global roughness 
exponent a distinct from the local one (the Hurst exponent H), previously reported [Ferreira et at, 
Appl. Phys. Lett. 88, 244103 (2006)]. A variety of scaling behaviors was obtained with varying 
substrate temperature. In particular, a transition from a intrinsically anomalous scaling regime with 
H 7^ a < 1 at low temperatures to a super-rough regime with H 7^ a > 1 at high temperatures 
was observed. The temperature is a growth parameter that controls both the interface roughness 
and dynamical scaling exponents. Nonlocal effects are pointed as the factors ruling the anomalous 
scaling behavior. 

PACS numbers: 68.55.-a,64.60.Ht,68.35.Ct,81.15.Aa 



I. INTRODUCTION 

The control of surface dynamics is a crucial step for 
production of thin film based optoelectronic devices since 
features as grain size and surface morphology are among 
the most important features affecting the efficiency of 
these devices.^ Consequently, the surface of many thin 
films has been studied extensively in the last years^H^ 
using dynamical scaling exponents which, in theoreti- 
cal studies, associate universality classes to distinct pro- 
cesses involved in the interface growthPSl In particular, 
the dynamical scaling of CdTe films, one of the most 
promising materials for the production o f high -efficiency 
solar-cells and other electronic devices ] 10 ^ 11 ! grown on 
amorphous substrates at varying temperatures has been 
characteri zed th rough the Hurst (H) and growth (f3) 
exponents! 4 * 7 * 12 ! CdTe films grown on amorphous sub- 
strates exhibited a peculiar behavior: the growth expo- 
nent and the global inte rfac e width are increasing func- 
tions of the temperature.™^ 

The analysis restricted to (3 and H exponents may con- 
tain only part of the information about the dy namical 
scaling, since several models^^ and experiment d 5 * 6 * 19 * 2 ^ 
have recently demonstrated the presence of anomalous 
scaling implying in a global roughness exponent distinct 
from the Hurst exponent, commonly called local rough- 
ness exponent. The interface evolution can be charac- 
terized by the interface width, defined as the root mean 
square deviation of the interface height around its mean 
value on a scale e defined by 



w(e,t) = i([h(x,t) - hf) 



1/2 



(1) 



where the bar represents the average inside windows of 
size e and (• • • ) the average over different profiles. The 
common behavior of the interface width of globally self- 
affine profile follows the Family- Vicsek scaling ansatz^^H 



w(e,t)=t?f 



(2) 



The scaling function f(u) is 
I const 



/(«) 



if u < 1 
if u > 1 



(3) 



where a is the roughness exponent. The horizontal corre- 
lation length grows with time as £ ~ t 1 ^ for sufficiently 
large substrates. In the literature, z is called dynamical 
exponent and is related to the growth and roughness ex- 
ponents bj0a = j3z. It is worth to stress that given the 
power laws w ~ and £ ~ t 1 ^ concomitantly with the 
scaling function |3]), the relation a = j3z must be obeyed 
and, therefore, these exponents are not independent. 

However, the scaling behavior of the local interface 
fluctuations may differ from the global ones character- 
izing an anomalous scalingP^H^ Indeed, the interface 
width on a scale e is given by 



w(e,t) 



-'Gb>) : 

where the anomalous scaling function is 

(u H if u < 1 
|_ const if u ^> 1 



Thus, the interface width is^Sl 



w(e,t) 



(4) 



(5) 



(6) 



where /3* = (a — H)/z. The standard self-affine Family- 
Vicsek scaling is recovered if a = H. 

In the present work, we report the dynamical scaling 
analysis of CdTe films grown on glass substrates cov- 
ered by fluorine doped tin oxide. In addition to the 
H and (3 exponents previously reported,™ we calculated 
the dynamic and roughness exponents, z and a. With 
the complete set of exponents it was possible to iden- 
tify new properties in this system, like anomalous scal- 
ing and super-roughness, which were not reported previ- 
ously. The paper is organized as follows. In section[TT|thc 
methodology is described. In section III the results are 
presented and discussed and, finally, some conclusions 
are drawn in section ITVl 
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II. METHODS 

t — 660 min 





200 



100 150 200 

Profile length (\im) 



(a) 





4111) 




200 




o 




4lll> 




200 











i 


4110 








200 


I 





u 




eh 


401) 




201) 

















400 




200 




1) 




4lll> 




21)1) 



_1_ 



50 100 150 200 250 300 

Profile length (|im) 

(b) " 

FIG. 1. Typical time sequence for scanned profiles of CdTe 
films grown at substrate temperatures (a) T = 150 °C and 
(b) T = 300 °C. 



CdTe films were grown on glass substrates covered by 
fluorine doped tin oxide by Hot Wall Epitaxy (HWE). 
CdTe films were produced with a growth rate of 1.4 A/s 
and growth times varying from t = 30 to 660 min at 
substrate temperatures varying from T = 150 to 300 °C. 
At least 20 surface profiles with length of 300 /im were 
measured for each sample, using a stylus profiler (XP1- 
AMBIOS) with vertical resolution better than 10 A and 
lateral resolution of 10 nm. Details of the experimental 
setup can be found elsewhere.^ 

In figure [IJ we show time sequences of scanned profiles 
of CdTe films illustrating the increasing of the interface 
width with substrate temperature previously reported.^ 



One can clearly beheld in figure [T] that growth instabili- 
ties characterized by sharp peaks emerge after a certain 
deposition time and are enhanced during the interface 
growth. Moreover, the instabilities are more pronounced 
at higher temperature. 

Growth and Hurst exponents for the samples used 
in this work were reported previously^ but the global 
roughness analysis was missing. Indeed, several exper- 
imental limitations prevent the production of samples 
with very long deposition times and, consequently, the 
usual relation w ~ L a for t — > oo cannot be usedP Alter- 
natively, the global roughness exponent a can be found 
applying the scaling relation a — f3z, since growth and 
dynamical exponents can be determined for short depo- 
sition times using the scaling laws w ~ t& and £ ~ t 1 ^ 2 , 
respectively. 

The correlation length was determined using the two- 
point correlation function T(e,i). For the sake of repro- 
ducibility, we used two distinct definitions to calculate 
r. The first one is the standard two point correlation 
function 



T(e,t) = (h(x + e,t)h(x,t) 



(7) 



where h is the detrended profile.^ The second one is 
the probability of the height difference between two sites 
separated by a distance e being lower than a fixed value 



r(e,i) = (Pr (\h(x + e,t)- h(x,t)\ < m 



(8) 



Where m is much shorter than the global interface width 
and much larger than the profile height resolution. The 
correlation function given by equation ([8]) is analogous to 
the two-particle correlation function commonly applied 
to fractal aggregates.^ Definitions (7| and ([8]) provide 
proportional correlation lengths and, consequently, the 
same dynamical exponent. Brackets and bars follow the 
same notation defined in equation ([!]). For the experi- 
mental data m = 0.1\h max — h m i n \ was used, where h max 
and h m i n are the maximum and the minimum heights in 
the detrended profile. 

In both definitions, T(e) was fitted by a two- 
exponential decay function 



ru.n : r„ + .1, cxp ( - M + A 2 ex P 



(9) 



where the fit parameters To, £i and £2 are time functions. 
The two-exponential decay is just a generalization of the 
standard exponential decay commonly used for correla- 
tion function fits. Our alternative choice is due to equa- 
tion (|9| fits very well all experimental data, while sin- 
gle exponential decay returns unsatisfactory fits for some 
samples. So, the correlation length is determined by 



I [T(M)-r ]de 
Jo 



= f / [r( e ) - r ]d£, 



(10) 
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where <C / < 1 is a fraction determining a characteris- 
tic decay. In the present work, we choose / — 0.1. 

In figure [2j the last method was applied to two distinct 
models in order to verify the validity of equation (pSl). In 
the first one, the Edwards- Wilkinson (EW) modelpf par- 
ticles are deposited on random sites of a linear chain and 
relax to the lowest neighboring position if this minimizes 
the particle height. In the second one, the Wolf- Villain 
(WV) modelp^ the particles are deposited at random 
and move to the nearest neighbor site that maximizes 
the number of bonds. As one can see, correct expo- 
nents were obtained for both models (see figure caption) . 
The method was also successfully applied to other well- 
established models. 




e (ML) 



FIG. 2. (Color on-line) Correlation function Q for the EW 
model at different deposition times. Symbols are simulational 
data and solid lines are non-linear fits using Eq. Q. Inset 
shows the correlation length as function of time for the EW 
and WV models. The slopes 1/z = 0.496 and l/z = 0.236 
for EW and WV models, respectively, are in very good agree- 
ment with the 1/2 and 1/4 expected valuesPflere, ML means 
monolayers. 



III. RESULTS 

Correlation functions using equation for T = 200° C 
at distinct times are shown in figure [3] and two represen- 
tative curves using definition (|8| are shown in figure[4j In 
both figures, the insets show the correlation length as a 
function of time for distinct temperatures, demonstrating 
the validity of the power law £ ~ t x / z . 

The dynamical exponent z as a function of tempera- 
ture is shown in the inset of figure [5j in which the error 
bars are the standard deviations of the exponents ob- 
tained for at least 4 groups of 5 profiles at each temper- 
ature and growth time. The exponent increases with the 
substrate temperature, implying in a slower spreading 
of the correlations along the surface and consequently an 
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FIG. 3. Correlation functions for the polycrystalline CdTe 
films grown at T = 200 °C using definition Q. Symbols are 
experimental data averaged over at least 20 surface profiles at 
each growth time and lines are nonlinear fits given by equation 
Q. In the inset, the time dependence of the correlation length 
is shown for two distinct temperatures. 
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FIG. 4. (Color on-line) Correlation functions for the polycrys- 
talline CdTe films using definition Q. Symbols are experi- 
mental data and lines nonlinear fits. In the inset, the time 
dependence of correlation length is shown for two distinct 
temperatures. 



increase of the interface width saturation time as temper- 
ature increases. The global roughness exponent was ob- 
tained from the scaling relation a = (3z using the values 
of the growth exponents (3 taken from a previous workP 
The results are shown in figure [5] Notice that the slow 
increase of the local roughness exponent H with tem- 
perature contrasts with the fast increase observed for a. 
Moreover, the system exhibits the so called intrinsically 
anomalous roughening,^ for which the local roughness 
exponent H < 1 is actually an independent exponent 
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and a may take values larger or smaller than 1. Indeed, 
theoretical arguments claim that symmetries and conser- 
vation laws restrict the emergence of intrinsically anoma- 
lous roughening in growth process controlled by nonlocal 
effects^ such as quenched disorder and shadowing ef- 
fects. It is very difficult to point out the nonlocal effects 
which could be ruling the growth dynamics of the sam- 
ples studied here, since they probably depend on details 
of the growth process. Theoretical studies and new ex- 
periments, including the investigation of other materials 
and substrates, are necessary to elucidate this point. 
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FIG. 5. Hurst and roughness exponents for CdTe films as 
functions of the substrate temperature. The inset shows the 
corresponding curves for the dynamical exponent. The values 
of the Hurst exponent were taken from Ferreira et aZP for 
t — 5h. Solid lines are nonlinear fits to guide the eyes. 

In figure [6j the lateral correlation length £y calculated 
using equation ^ is shown as a function of temperature 
for a fixed time. For low temperatures, the correlation 
length £j| increases with temperature as it is expected 
since the larger the diffusion rates the larger the corre- 
lation spreading. However, the opposite is observed at 
higher temperatures, which reflects the onset of peaks, 
i.e., the transition to the regime ruled by instabilities. 
The microscopic origin of this behavior is not clear. Ex- 
ponents obtained from correlation functions Q and ([8| 
diverge at less that 10 % for any temperature, assuring 
the reproducibility of the results. 

IV. CONCLUSIONS 

CdTe films grown on glass substrates covered by fluo- 
rine doped tin oxide by Hot Wall Epitaxy (HWE) were 
studied using the interface dynamical scaling theory. We 
determined the dynamical exponent which revealed an 
intrinsically anomalous scaling characterized by a global 
roughness exponent a distinct from the local one, the 
Hurst exponent H. In addition to the control of surface 
width formerly reported,^ we showed that the scaling ex- 
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FIG. 6. Correlation length as a function of temperature for 
distinct times. 



ponents and, consequently, the growth dynamics in the 
surface, can be modified by varying the control parame- 
ters, in particular the substrate temperature. The tran- 
sition from a < 1 to a > 1 has a particular importance 
which can be observed in figure [TJ For T — 150° C, when 
a < 1, the interface fluctuations grow approximately uni- 
formly along the surface, while for T — 300°C, when 
a > 1, the surface is ruled by instabilities. 

As discussed by Lopez et aij^ symmetries and con- 
servation laws restrict the emergence of the intrinsically 
anomalous roughening in local growth models. Indeed, 
disorder and/or nonlocal effects are required for intrin- 
sic anomalous roughening. Our experiments are consis- 
tent with this hypothesis if we suppose that the original 
substrate introduces a sort of nonlocal disorder, which 
amplifies the instabilities along the surface growth. It 
was observed in simulations of epitaxial growth with a 
quenched disorder in the step barriersj^ that the growth 
exponent /3 is an increasing function of temperature. In 
addition, the shadowing effects are probably strong in 
HWE technique due to the collimation vapor directed 
to the filmP In summary, our results are in accordance 
with the conjecture that the emergence of the intrinsi- 
cally anomalous roughening in local growth process re- 
quires the break of symmetries and/or conservation laws. 
Following the simulational results by Elsholz et alP, we 
propose that the nonuniformity of the original substrate 
introduces a quenched disorder as the cause of symme- 
try break, which is enhanced by a nonlocal shadowing 
effects came from the HWE deposition method. Theo- 
retical studies and new experiments using other materials 
and substrates are necessary to elucidate this point. 
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